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Abstract 

It is proved in this paper that Riemann-Liouville processes can arise from the 
temporal structures of the scaled occupation time fluctuation limits of the site- 
dependent (d, a, cr(x)) branching particle systems in the case ofl=e?<a:<2 and 
J R a(x)dx < oo. 
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1. Introduction 

A Riemann-Liouville process (RL process) is a centered Gaussian process. It can be 
represented by the moving average representation of Brownian motion 

t 



(t-u) H ~ 1/2 B(du), 0<t<oo. (1.1) 



o 



Here B(t) denotes the standard Brownian motion and H > is a constant. For 
convenience, we denote the RL process by {R H (t),t > 0}. It is easy to see from 
(jl.ip that R H does not have stationary increments but is self-similar with index H. 
RL processes are very useful and important because they are closely related to the 
Riemann-Liouville derivative of Gaussian noise. As early as 1950s, Levy [11] had briefly 
commented on this kind of processes. Since then, much work has been done in order 
to study it; see, for example, []Jll [15j [18] and the references therein. 

From (|l.ip one can readily think of the fractional Brownian motion (FBM) B H with 
(Hurst) index H G (0, 1). The FBM B H has a moving average representation as follows. 

B H (t) := f \(t - n)J" 1/2 - (-n)J" 1/2 j B(du), (1.2) 

where u+ = u V 0. It is well-known that B is a centered Gaussian process with the 
covariance function 

E[B H (t)B H (s)} = ±(\s\ 2H + \t\ 2H -\t- s\ 2H ). (1.3) 
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and B H is self-similar with index H and has stationary increments. 

Sub-fractional Brownian motions (Sub-FBMs) were formally introduced in 2004 in- 
dependently by Bojdecki et al [3] and Dzhaparidze and van Zanten [9], respectively. 
The Sub-FBM {rj H (t),t > 0} with index H G (0, 1) is also a centered Gaussian process 
with the covariance function 

E[ V H (t) V H (s)] = \s\ 2H + \t\ 2H - + t) 2H + \t- s\ 2H ). (1.4) 

This kind of processes has many properties which are the same as those of FBMs, but 
its increments are not stationary. 

When H ^ 1/2, RL processes, FBMs and Sub-FBMs share the same property that 
their increments on non-overlapping intervals are correlated. In fact, Bojdecki et al [3] 
showed that the covariance of increments on intervals at distance t decays like t 2H ~ 3 
for the Sub-FBM rj H and t 2H ~ 2 for the FBM B H . In addition, as we will point out in 
Proposition 2.1, this index for the RL process 

R H ig t H-3/2 _ 

Bojdecki et al O H] studied the functional limits of the occupation time fluctuation 
processes 

1 f Tt 

X T (t) = — / (N(s) - E(N(s)))ds, t > 0, (1.5) 
Ft Jo 

of a classical (d, a, l)-branching particle system, where N(t) denotes its random count- 
ing measure at time t, and found that under the assumptions a < d < 2a and 
Ft = Tz~2E, a Sub-FBM rj H with H = | — can arise from the temporal structure 
of the limit process of Xr(t). Here E(iV(s)) is the expectation function understood as 
(E(N(s)), <f>) = E((N(s), <j>)) for any <f> G S{R d ), the space of smooth rapidly decreasing 
functions. Motivated by the aforementioned works, the main purpose of this paper is to 
report a similar but new result. We show that under certain conditions, RL processes 
arise from the limits of occupation time fluctuations of the so-called site-dependent 
branching particle systems (see Li p3] for details). This result is not obvious at a first 
sight. In fact, so far as we know, it is not formally reported in existing literature even 
if RL processes have been explored for longer time than that of FBMs and Sub-FBMs. 

Precisely, in this paper, we consider a branching particle system in M, whose particles 
start off at time t = from a Poisson random field with Lebesgue intensity measure A, 
evolve independently with a symmetric a-stable Levy motion and undergo at rate 7 the 
critical, finite variance branching which depends on the position x where the particle 
splits and is controlled by the moment generating function 

g(s,x) =s + <t(x)(1 - s) 2 , < a < 1,0 < a(x) < 1/2, x G M, 

where cr(-) is a measurable function on R. For convenience, we refer this model to 
the (1, a, cr(x))-particle system. Let N(s) denote the random counting measure of the 
(1, a, cr(x))-particle system. Consider the scaled occupation time fluctuation process 
defined by (jl.5p . We prove that when a G (1,2), and Ft = T5"o i the RL process R H 
with H = I — — arises from the temporal structure of the fluctuation limit process if 
f^a(x)dx < 00. In addition, we also prove that the standard Brownian motion can 
arise from the limit process in the case of a = 1 if Ft = T 1 / 2 InT and L a(x)dx < 00. 
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There is much other literature in the field of functional limit of occupation time 
fluctuations of branching systems. For example, Li et al [12|, [T3] studied the branching 
particle systems with degenerate branching. Milos [16] investigated the change caused 
by the varied initial distributions of the classical (d, a, j3) branching particle systems. 

Without other statement, in this paper, we use M to denote an unspecified positive 
finite constant which may not necessarily be the same in each occurrence, and use 
B, B H ,n H , R H to denote the standard Brownian motion, the FBM with index H, the 
Sub-FBM with index H and the RL process with index H, respectively. 

The remainder of this paper is organized as follows. In Section 2, we report the main 
results. Section 3 contains some basic calculations. In the last section, i.e. Section 4, 
we prove the main results. 

2. Main results 

From now on, we always let N(s) and -Xr(t) be the random counting measure of the 
(1, a, cr(a;))-particle system and the corresponding scaled occupation time fluctuation 
process. Let £ denote the particles' space motion. Then £ is the symmetric a-stable 
Levy process. Denote its semigroup by {L t }t>o and the transition density by pt, i.e., 

L t f{x) :=E(/(£(t + a))|£( S ) = x)= / p t (y - x)f(y)dy, 

JR 

for all s, t > 0, x G R and bounded measurable functions / (to avoid ambiguity we 
sometimes write L t f(x) as L t (f(-))(x)). For every (p £ <S(R), denote the Green potential 
operator by 

;>oo 

Gtp(x) = / L s ip(x)ds, 
Jo 



and define 



From Li [H] , 



G t ip{x) = I L s ip(x)ds. 



E((JV(t),0))= / L t( f ) {x)dx = / <j>{x)dx = (\4>). (2.1) 
Jr Jr 

The occupation time fluctuation process Xt = {Xx{t),t > 0} is rewritten as follows. 



1 



Tt 



(X T (t), <P) = T r (N(s) - A, fids, (2.2) 

b T Jo 

for every <j) S <S(M), where Ft is a suitable norming constant. 
The main results of this paper read as follows. 

Theorem 2.1 Suppose D := j R a(x)dx < oo and a £ (1,2). Let F T = T 3 ~ 2 / Q . Then 
X T => K\R H in C([0,t],R) for any r > 0, where R H {-) is a RL process with index 
H = 3/2 - l/a and K = ^DT{1/ 'a) /[vr(a - 1)]. 

Theorem 2.2 Suppose D := f R a(x)dx < oo and a = 1. Let F T = T(lnT) 2 . 
Then Xt =>i CXB, where C = Here means that for each t > 0, 

Jq (Xt, s))ds converge to (X, ^(-, s))ds in distribution for all if) € 5(M 2 ). 
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The long range dependence of processes can be characterized by means of the so-called 
dependence exponent which was developed in [6]. According to [6], the dependence 
exponent k of a Gaussian process {Ct]t > 0} is 

k = inf sup I 7 > : Cov(C„ - ( u , C, T +t - Ct+s) = o(T~ 7 ) as T ->■ oo \. (2.3) 

0<u<v<s<t K ) 

The following proposition is of the dependence exponent of RL processes. It seems 
simple; we write down it because we do not find an appropriate reference. 

Proposition 2.1 The dependence exponent k of the RL process {R H (t),t > 0} with 
\ %s\-H. 

Below we make some comments on our results. 

Remark 2.1 (1) Prom the assumption of J R a(x)dx < oo, one can readily relate the 
(1, a, u(x))-particle systems to the particle systems without branching. Some results 
of the latter 's occupation time fluctuations were revealed in Bojdecki et al [H[5]. They 
considered the particle systems where particles start off at time t = from a Poisson 
random field with Lebesgue intensity measure A and move independently in R d with 
a symmetrical a-stable Levy motion, and proved that if d/a < 1, then Ft = T 1- ^ 
and the limit process X = KXB 1 ^^; if d/a = 1, then Ft = (TlnT) 1 / 2 and the limit 
process X = KXB. Comparing their results with the corresponding results in this 
paper, there at least exist two differences in the case of d = 1. They are, respectively, 
F T = T 1 / 2 InT vs F T = (TlnT) 1 / 2 when a = 1, and F T = ri~« plus R^~* vs 
Ft = T 1_ 2^ plus B 1 "^ when a & (1, 2). The limit processes Rz~a and B 1 "^ 
have the same dependence exponent 1/a but different self-similarity indice: | — — vs 
1 — 2o> an d different increments' property: non-stationary increments vs stationary 
increments. The difference on the norming Ft also shows that the occupation time 
fluctuations in the (1, a, cr(x)) model increase faster than those in the particle systems 
without branching. 

(2) It is also interesting to compare the results of (1, a, <j(x))-particle systems where 
o~(x) = 1/2 with the results in this paper. In fact, when cr{x) = 1/2, the (1, a, o~(x))- 
particle systems are the classical (d, a, /3)-branching particle systems with d = 1, (3 = 1 
(see, for example, Bojdecki et al [4]). For the classical (1, a, l)-branching particle 
systems, Theorem 2.2 in Bojdecki et al [1] (see, also, Bojdecki et al ) told us that a 

3 1 

long-range dependent Sub-FBM rj2~2E with dependence exponent 1/a arise from the 
temporal structure of the limit process of occupation time fluctuations in the case of 
a G (1/2, 1). Moreover, Bojdecki et al intuitively explained the cause of the long range 
dependence as the "clan recurrence". In this paper, we get the long-range dependent 
RL process R^~a in the case of a £ (1)2). The long-range dependence exponent is 
1/a < 1, and the long-range dependence is likely caused by the recurrence of particles' 
motion. 

(3) By the terminology in Bojdecki et al [7], the convergence in Theorem 2.2 
is called convergence in the space-time, or integral, sense. Bojdecki et al [2] pointed 
out that convergence resembles the convergence of finite-dimensional distributions 
but neither implies the other. However, by the same proof as in Bojdecki et al [7j, we 
can additionally conclude in Theorem 2.2 that Xt converges to CXB in the sense of 
finite-dimensional distributions. 
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At the end of this section, let us show some heuristics to explain why one can expect 
the RL process arising as a limit in Theorem 12.11 To this end, observe the covariance 
of Xt- Let < s, t < 1. From (|4,2p . it is easy to get that 

Cov((X T (t)^),(X T (s),^)) = ^ f f Cov({N(Tu)A),{N{Tv)A))&udv 



J JO 



<t(x) / L r ( M _ r )<^(x)L< r („_ r )<^(x)drdxdwi'y 



JO JO 

^2 rs rt 



+^2 f f f (p(x)L Tlu _ vl< j)(x)dxdudv. 
Jo Jo Jr 



T JO JO 

From (3.1) and (3.2) in Bojdecki et al [6], for any <p G <S(M), there is constant M > 
Mi > such that 

Mru- 1 ^ < L u (j)(x) < Mu~ 1/a (2.4) 

for any u > 0. Since = T 3 ~q and 1 > 1/a, we have that 
r s ft r ruAv 



F 2 

T 



rs rt p pu/\v 

/ / / a(x) / L T ^ u _ r - ) if>(x)L T ^ v _ r - ) (f)(x)drdxdudv 
Jo Jo JR Jo 

ruAv 

dudv / {u-r)- l/a {v - r)~ 1/a dr 
Jo 

( s _ r )i-V«( t _ r )i-i/« dr) 



JO 

sAt 



where A ~ B means that A/B S (0,oo), and that 

-Q. rs ft 



rpz rs rt r 

lim —2 / / / 4>(x)L T i u _ v i(f>(x)dxdudv = 0. 

Therefore, 

rsAt 

Cov((X T (t),^), (X T (s), <P)) ~ / (a - r) 1 " 1 /^ _ rf^r, 

Jo 

which is the covariance of the RL process R H with H = 3/2 — 1/a. 

3. Preliminary calculations 

Define a sequence of random variables Xt in 5'(R 2 ) as follows: For any T > 1 and 
V> E 5(M 2 ), let 

(X T ,V) = f\x T {t),^{;t))dt. (3.1) 







Without other statement, in the sequel, ip € 5(R 2 ) always has the form ip(x,t) 
(f)(x)h(t), where <fr E <S(1R) and /i E S(R) are nonnegative functions. Let 



1 1 ,, .?,s 



h(s)= I h(t)dt and ih(x,s) = — Hx)h(-). (3.2) 
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Define 



V^ T {x,t,r) = 1-E x (exp{ -J {N{s), ip T (-, r + a))da}) . 



(3.3) 



where K x (f(N)) := ¥,(f(N)\N(0) = e x ) and e x denotes the unit measure concentrated 
at x G M. From Section 3 in p3], we get that 



V^ T (x,t,r) = / L s [if) T (-,r + s)(l - V^, T (-,t - s,r + s))](x)d& 
Jo 



-7 / L s [a(-)VL(;t-s,r + s)](x)ds. 



Furthermore, define 



Then 



and 



where 



and 



E(e 



-{x 7 



h(T. 



J^ T (x,t,r) := L s ifj T (-,r + s)(x)ds. 
Jo 



Vip T (x,t,r) < Jip T (x,t,r), 
-- exp (h(T, fa) + I 2 (T, fa) + I 3 (T, 

T — s, s)ds, 



\ih) = 7 / cr(x)dx / J% T (x, 
JR JO 



/ 2 (T,Vt) = 7 / [ T [Vj T (x,T 

JR Jo 



s,s) - Jl Jx,T - s,s)]ds 



> -27 



hi(T, ip T ) + 7/22 (T 1 , V't' 



I 3 (T,if) T )= dx ip T (x,s)V^ T (x,T - s,s)ds, 
Jr Jo 

hi(T,ipT)= / cr(x)dx / J^ T (x,T — s, s)ds 
Jr Jo 

f T 

x / L u ^ s (Tp T (-,u)J^ T (-,T - u,u))(x)du, 

J s 

l22(T,ifar)= / cr(x)dx / Jip T (x,T — s,s)ds 
Jr Jo 

x / L u - S (a(-)J^ T (-,T - u,u))(x)du. 



(3.4) 

(3.5) 

(3.6) 
(3.7) 

(3.8) 

(3.9) 
(3.10) 



(3.11) 



(3.12) 



3.1 Limits of I\, I2 and ^3 when 1 < a < 2 
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Step 1. We consider the limit of I\{T,iPt)- From (|3.8p and (|3.5p . we get that 
h(T, ^ T )=7 / a(x)dx / I / L u tp T (x, s + u)du I ds 

JO \J0 J 
pi pi prAt 

= / fc(t)dt / /i(r)dr / * T (fl,r,t)dfl, (3.13) 
io Jo jo 



where 



T(t-s) pT(r-s) 



^ T (s,r,t) = — 1 I a(x)dx I L u (f)(x)du I L v qb{x)dv. (3.14) 

Ft Jr jo jo 

The self-similarity of the a-stable Levy process £ implies that 

-ir- r L u f{x)du= f f u- 1 / a p 1 ((x-y)T- 1 / a u- 1 / a )f(y)dy, 
T a Jo Jo Jr 

for any integrable function /. This means that for every < t < 1 and non- negative 
integrable function / 

a i 



T i 1_ «Pi(0) / f(y)dy> — - / L u /(x)dx. (3.15) 
1 Jr T « Jo 



a 

Additionally, as T — > oo, 
T 

Therefore, (|3.14p and the facts that J K a(x)dx < oo and F T = T 6 ~ 2 / a lead to 

2 



-Lr L u f(x)dx^^—t l -^ Pl (0) [ f(y)dy. 
Jo a ~ 1 Jr 



9 T (a,r,t)^(p 1 {0)-^) 2 [(t-s){r-8)] 1 - 1 / a / a{x)dx( / 0(s)d*Y 
a - 1 Jr v Jr 

= C[(t-s)(r-s)] 1 - l / a (f ^(x)dx)\ 



as T — t- oo, where 

/or(o)\2 /■ / r(i/a) \ 2 /■ , v , 

C = 7( / a{x)dx = j[— -I / a(x)dx. 

\ a - 1 / J R Vvr(a -I) J J R 

Consequently, from (|3. 13|) we can readily have that as T — > oo, 

h(T,ip T ) ^C( [ (f>(x)dx) 2 [ h{t)dt [ h(r)dr [ [(t - s)(r - s^^ds. (3.16) 
Wr 'Jo Jo Jo 

Step 2. We consider the limit of Ii(T,i\)t)- Using (|3.2|) . (|3.5|) and the fact that /i is 
bounded, we have that 

M 

Jjj T (x,T - u,u) < — G T 4>(x), (3.17) 
ft 



Yuqiang LI 



for any u < T and some M > 0. Applying this inequality to (|3.1ip and using (|3.2p 
again, we get a constant M > such that 

MT /" /" T /" T /" T 

hi(T,ipr) < / cr(x) / L v (p(x)dv / L u (0 / L w (j)dw)(x)dudx 

JO JO JO 

MT f 

<— 3- / a{x)G T (t){x)GT{(t)GT^){x)dx. (3.18) 
Ft j^ 

It follows from ()3. 15[) that for any non-negative integrable function / there exists a 
constant M > such that 

G T f(x) < MT l ~ l l a . (3.19) 
Therefore, substituting F| = T 3 ~ 2 / a into (IBTTHD . we can readily obtain that 

M f 

hi (T, ifr) < a(x)dx -> 0. (3.20) 

Furthermore, applying (|3.17p and (I3.2p to (|3.12p . we obtain that 

MT f 

/22(T,*/T)<-=r / ^(^)(GT^)) 2 GT(cTG T ^(x)dx, (3.21) 

for some M > 0. From (|3.19p . (|3.2ip and the fact J R a(x)dx < 00, it follows that 

M 

h2{T^ T )< Tl/a _ 1/2 ->0, (3.22) 

where the convergence is due to the fact a G (1,2). Note that hiT^tpj') < which 
follows from (|31)]) . Combining (|3U]) with (f3TT8|) and (f3722|) yields that as T -)• 00, 

/ 2 (T,^r)->0. (3.23) 

Step 3. At last we pass to the limit of I$(T,ipT)- Let 

hi{T,ipr) := dx ^t(x, s)J^ t (x, T - s, s)ds. (3.24) 

j 

From (fBTTUD . (fBTol) . (fHTTTD and (pD^D . it follows that 

MT f 

< h{T,i> T )<hx{T^ T ) < -=g- / <l>(x)G T <l>(x)dx 

< =g / ^d*. ( 3 - 25 ) 

P T JR 

for some M > 0. Substituting F| = T 3 ~ 2 / a into (|335jl . q G (1, 2) indicates that 

/ 3 (T,^ T )^0. (3.26) 
3.2 Limits of I\, I2 and I3 when a = 1 
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We first point out that 



and 



lim h{T^ T ) = 0, (3.27) 

T->oo 



lim I 3 (T,ipr) = 0. (3.28) 

T->oo 



The details are similar to those of I2 and I3 in Section 3.1, and hence we omit them. 
Below, we discuss Ii(T,iJ>t)- 

From (3.47) in [8] we know that for any < t < 1, 



1 



Tt 



lim — / L u cj>(x)du = pi(0) / <f>(x)dx. (3.29) 

Using this fact and noting that J K a(x)dx < 00 and Fj, = T(lnT) 2 , we get that 

lim ^t(s, r, t) = 7 pi(0) f (p(x)dx j a(x)dx, (3.30) 

where ^x(s,r,t) is the same as (|3.14p . Note that (|3. 19[) and the integrability of a(x) 
imply the boundedness of *$>t(s, r, t) with respect to T > 0, < r, t < 1 and < s < rAt. 
Therefore, the dominated convergence theorem and ()3.13|) imply that as T — > 00, 



where 



h(T,i/; T ) -»■ C7i( / (f)(x)dx) / h(t)dt h{r){r M)dr, (3.31) 

J M J •/ 

d = 7(Pi(0)) 2 [ (J(x)dx = -1 / a(x)dx. 
4. Proofs of the main results 

Proof of Theorem 2.1 Without loss of generality, fix r = 1. As explained in Bojdecki 
et al [l],to prove this theorem it is sufficient to show that 

(i) For each non-neg ative tp G <S(M d+1 ), as T — > 00 

E(e-<*^>) -> exp Q jf jf Cov((X( S )>(-, a)), <X(t), </>(", t)))dsdt) ,(4.1) 

where X n and X are defined as (I3.ip and X is the corresponding limit process; 

(ii) {{Xt, (f));T > 1} is tight in C([0, 1],K) for any given (j) G <S(K), where the theorem 
of Mitoma [17] is used. 



If ip(x,t) = (p(x)h(t), where (j) G S(M) and h G <S(M) are nonnegative functions, then 
the limit process in Theorem 2.1 satisfies that 

Cov((X( S ), ^(-, 8 )>, (X(t), </>(•, t)>) = K 2 (A, 0) 2 /i(8)/i(t)E(l?(a)12(t)), 

/•sAt 

= K 2 (A,^) 2 / i (s)/ l (t) / [(t-ujfs-u)] 1 - 1 /^!*, 
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where K 2 = 27 J R a (x) dx ( j^rr^ ) 2 ■ Furthermore, combining (|3.7p with (|3.23[) . ()3.16|) 
and (^261) yields 



lim E(e- (Xr 'W) 

T— »oo 



exp 7 



(i C^^f} ) 2 I v(x)dx C f h(s)h(t) [(t - u)(s - u^-^dudsdt] 
V V vr(a - 1) / J R Jo J J ) 

= expQjT y Cov(<X( S ),V(-,s)),<X(t),^(-,t)))d S diY 

This means that (|4.ip holds for the special case ip(x,t) = 4>{x)h{t). 

For general ip, the proof of (i) is the same with slightly more complicated notation. 
The details are omitted. 

We now pass to prove (ii). Note that by some simple and standard arguments, one 
can readily get that for any f,g& S(M) and s <t, 

E x ((N(s),f)(N(t),g)^ = L s (fL t _ s g)(x) + 2 7 J\ u (aL s _ u fL t _ u g)(x)du. (4.2) 

Thus, by Poisson initial condition, 

Cov({N(s)J),{N(t),g)) = [ f(x)L t - s g{x)dx+2 1 I a(x) [ L s ^ u f(x)L t _ u g(x)dudx. 
v ji Jr Jo 

From the setting of Xxit), we have that for any < s < t < 1, 

E((X T (t),ct)) - {X T (s)^)f = ^2^ f s Cov^(N(Tu),(j>),(N(Tv)^))dudv 
27T 3 rt r ruAv 



p2 

rr p Js Js 
^2 r t rt 



a(x) / L T ( u _ r -)4>(x)L T ( v _ r }(/)(x)drdxdudv 



+ [ [ [ 4>{x)L T \ u _ v \4>(x)dxdudv. (4.3) 

*T J s J s Jr 

Applying the last inequality of (|2.4p to the last two formulas of (|4.3p . noting that 



= T 3 a and cr{x) is integrable, we get that for all T > 1 
27T 3 /"* /"* f f uAv 

—EvT III a ( x ) / L T(u~r) <t>(x)L T (v-r) <t>{x)drdxdudv 
"t Js Js Jr Jo 

pt ft puAv 

<M du dv {u-r)- l / a {v -r)~ l / a dr 

Js Js Jo 

<M [ du I \u-v\~ 1/a dv 

-1 rl 

'0 Jo 



M(t-sf-~ f [ \u-v\~ l / a dudv < M{t-s) 2 ~«, (4.4) 
Jo Jo 



and 



2 H rt 



-p2 

r T J s J s 



4>{x)L T \ [U _ v \ [ (f){x)dxdudv 



< ' - / \u-v\- l/a dudv<M(t-s) 2 ~^, (4.5) 



10 



Riemann- Liouville processes and Branching Systems 



11 



hold for constants M > which only depend on (p. 
Combining (|Pj) with (fO)) and (f43j) . we get that 

E«Xr(t), 0) - <A T (s), 0)) 2 < M(0)(t - s 



Therefore, from Billingsley [U Theorem 12.3] and the fact (Xr(0),^>) = 0, it follows 
that {(Xt, 4>)}t>i is tight in C([0, 1],M). The proof of Theorem 2.1 is complete. □ 



Proof of Theorem 2.2. The proof of Theorem 2.2 is easy. In fact, combining ([37 
with (|3.27|) . (|3.28|) and (|3.31|) we can readily get the corresponding formula (|4.ip . More 
details are omitted. □ 

At last, we give the proof of Proposition 2.1. 

Proof of Proposition 2.1. From the representation (jl.ip of R H (■), it is easy to see 
that for any 0<u<v<s<t, 

Cov(R(v) - R(u), R(T + t) - R(T + s)) 

/•OD 

= / [(«-r) ff -a {r < w} -(u-r) ff -U {r < u} ][(r + t-r) H -s-(T + a -r)^]dr. 
J o 

By Taylor expansion formula, it is easy to see that for any given t,s,r, 

• p., t , ., r) [(r + t-r)»-i -(T + .-r)-i] ^ _ 1 _ 

T 2 ^ 

as T — > oo. Consequently, one can readily have that 

Cov(i?(i;) - R(u),R(T + t) - R(T + s)) 



lim 



1 f°° 1 1 

= (H- -){t - s) [{v - r)^-3l {r < w} - (u- r) H -2l {r < u} )&r 

= 2 ^\(v H ^-u^)(t-s), (4.6) 

which and (|2.3p together imply that the dependence exponent of RL process 

with H^hial-H. □ 



Acknowledgements 

The author thanks the anonymous referees for careful reading of the paper and for 
the useful suggestions which have helped to improve the paper significantly. 



References 

[1] Billingsley, P. Convergence of Probability Measures, Wiley, New York, 1968. 

[2] Bojdecki, T., Gorostiza, L, and Ramaswami, S. Convergence of iS'-valued processes and 
space-time random fields. J. Fund. Anal., 66 (1986): 21-41. 

[3] Bojdecki, T., Gorostiza, L, and Talarczyk, A. Sub-fractional Brownian motion and its 
relation to occupation times. Statist. Probab. Lett., 69 (2004): 405-419. 



11 



12 Yuqiang LI 



Bojdecki, T., Gorostiza, L, and Talarczyk, A. Limit theorem for occupation time fluc- 
tuations of branching systems I: Long-range dependence. Stochastic Process. Appl., 116 
(2006): 1-18. 

Bojdecki, T., Gorostiza, L, and Talarczyk, A. Limit theorem for occupation time fluctu- 
ations of branching systems II: Critical and large dimensions. Stochastic Process. Appl., 
116 (2006): 19-35. 

Bojdecki, T., Gorostiza, L, and Talarczyk, A. A long range dependence stable process and 
an infinite variance branching system. Ann. Probab., 35 (2007): 500-527. 

Bo jdecki, T., Gorostiza, L, and Talarczyk, A. Occupation time fluctuations of an infinite- 
variance branching system in large dimension. Bernoulli, 13 (2007): 20-39. 

Bojdecki, T., Gorostiza, L, and Talarczyk, A. Occupation time limits of inhomogeneous 
Poisson systems of independent particles. Stochastic Process. Appi, 118 (2008): 28-52. 

Dzhaparidze, K.O. and van Zanten, J.H. A series expansion of fractional Brownian motion. 
Probab. Theory Related Fields 130 (2004): 39-55. 

Hclga, S. An optimal wavelet series expansion of the Riemann-Liouville process. J. Theoret. 
Probab., 22 (2009): 1030-1057. 

Levy, P. Random functions: General theory with special reference to Laplacian random 
functions. Univ. California Publ. Statist., 1(1953) 331-390 

Li, Y. and Xiao, Y. Occupation time fluctuations of weakly degenerate branching systems. 
J. Theoret. Probab. DOI: 10.1007/sl0959-011-0358-3 

Li, Y. Occupation time fluctuations of strongly degenerate branching systems. Sci. China 
Ser. A Mathematics. 54 (2011): 2669-2682. 

Li, Y. Fluctuation limits of site-dependent Branching systems in critical and large dimen- 
sions. Stat. Probab. Lett., 81 (2011): 1604-1611. 

Lifshits, M. Linde, W. and Shi, Z. Small deviations of Riemann-Liouville processes in 
L g -space with respect to fractal measures. Proc. London Math. Soc, 92 (2006): 224-250. 

Milos, P. Occupation time fluctuations of Poisson and equilibrium finite variance branching 
systems. Probab. Math. Statist., 27 (2007): 181-203. 

Mitoma, I. Tightness of probability on C([0, l\,y') and D([Q, Ann. Probab. 11 

(1983), 989-999. 

Sithi, V.M. and Lim, S.C. On the spectra of Riemann-Liouville fractional Brownian mo- 
tion. J. phys. A Math. Gen., 28 (1995): 2995-3003. 



12 



